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The electron motion along a chain is described by a continuum version of the Su-Schrieffer-Heeger 
Hamiltonian in which phonon fields and electronic coordinates are mapped onto the time scale. The 
path integral formalism allows us to derive the non local source action for the particle interacting 
with the oscillators bath. The method can be applied for any value of the e-ph coupling. The path 
integral dependence on the model parameters has been analysed by computing the partition function 
60! and some thermodynamical properties from T = IK up to room temperature. A peculiar upturn 

in the low temperature heat capacity over temperature ratio (pointing to a glassy like behavior) has 
been ascribed to the time dependent electronic hopping along the chain. 
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C/3 1 I. INTRODUCTION 

o5 . 

Conjugated organic polymers have been a focus of 
G ■ intensive research over the last twentyfive years after 
' the discovery of their remarkable conducting properties 
tuned by the dopant concentration [1]. Undoped trans- 
polyacetylene sustains neutral solitons in the dimerized 
O , chain as a consequence of the Peierls instability [2]. 

Charge injection in polymers induces a lattice distortion 
t—{ | with the associated formation of localized excitations, po- 
J> ■ larons and/or charged solitons. While the latter can exist 
only in iraras-polyacetylene whose ground state is twofold 
degenerate, the former solutions are more general since 
they do not require such degeneracy. The Su-Schrieffer- 
Heeger (SSH) model Hamiltonian [3] is the fundamental 
tool in the analysis of quasi one dimensional polymers. 
In the weak coupling regime a continuum version [4] of 
■ the SSH model has been developed and polaronic solu- 
tions have been obtained analytically [5]. Although the 
ID properties of the SSH model have been mainly in- 
vestigated so far [6-11], extensions to two dimensions 
were considered in the late eighties in connection with 
the Fermi surface nesting effect on quasi 2D high T c su- 
perconductors [12]. Recent numerical analysis [13] have 
revealed the rich physics of 2D-SSH polarons whose mass 
seems to be larger than in the ID case, at least in the 
intermediate regime of the adiabatic parameter [14]. 

One of the main issues in the study of e-ph models 
regards the possibility to monitor the conditions for po- 
laron formation as a function of the coupling strength. 
The Feynman path integral method [15] is a landmark in 
this respect and several studies [16-18] have applied the 
path integral formalism both to the acoustical polaron of 
the Frohlich Hamiltonian and to the Holstein model [19]. 
As a main property of the SSH model, the e-ph inter- 
action modifies the bare electron hopping integral of the 
tight binding approximation thus leading to an Hamil- 
tonian linearly dependent on the phonon displacement 
field. This feature will allow us to attack the SSH model 
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by path integrals techniques after introducing a general- 
ized version of the semiclassical model which treats only 
the electrons quantum mechanically. The general formal- 
ism is described in Section II together with the method 
to compute the full partition function of the interacting 
system. In Section III, we focus on the peculiarities of the 
free energy and heat capacity as a function of the model 
parameters. Section IV contains the final remarks. 

II. PATH INTEGRAL FORMALISM 

The ID SSH interacting Hamiltonian reads: 

H = ^ Jr,r+l(/ r f /r+l + fr+lfr) 
r 

J r ,r+1 = -\[ J + a ( U r - Ur+l)] (0.1) 

where J is the nearest neighbors hopping integral for 
an undistorted chain, a is the electron-phonon coupling, 
u r is the dimerization coordinate which specifies the dis- 
placement of the monomer group on the r— lattice site 
along the molecular axis, /J and f r create and destroy 
electrons (i.e., tt band electrons in polyacetylene) on the 
r— group. Our non-interacting Hamiltonian is given by a 
set of independent oscillators which will be treated clas- 
sically in the following. The real space Hamiltonian in 
cq.(l) can be transformed in a time dependent Hamil- 
tonian by introducing x(t) and y(r') as the electron co- 
ordinates at the r and r + 1 lattice sites, respectively. 
The spatial e-ph correlations contained in eq.(l) are then 
mapped onto the time axis by changing: u r — » u(t) and 
u r +i — > u(t'). Accordingly, the time dependent Hamil- 
tonian is: 

H(t,t>) = J T}T ^p(x(r))f(y(r')) + p(y(T'))f(x(r))) 
Jt,t> = -\[J + a(u(r) - u(t'))] (0.2) 

The twofold degenerate ground state of the SSH Hamil- 
tonian undergoes a Peierls instability which opens a gap 
at the Fermi energy. In real space the soliton connects 
the two degenerate phases with different senses of dimer- 
izations and a localized electronic state is associated with 
each soliton. Both electron hopping between solitons 
and electron hopping to band states (thermal excita- 
tion) are in principle possible within the model. Mapping 
the Hamiltonian onto the time scale we maintain these 
properties while thermally activated electron hops be- 
come time dependent. H (r, r') is therefore more general 
than the real space SSH Hamiltonian since hopping pro- 
cesses are not constrained to first neighbors sites along 
the chain. We emphasize that a variable range hopping 
introduces some local disorder in the system. 

Eq.(2) displays the semiclassical nature of the model 
in which quantum mechanical degrees of freedom interact 
with the classical variables u(t). Setting r' = 0, u(0) = 
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?/(0) = and averaging the electron operators over the 
ground state we obtain the time dependent semiclassical 
energy which is linear in the displacements: 

<H(t) >= -^(j + au(r))(G[x(T),T}+G[-x(T),-T]) 

L r /a dk ^exp(-iv n T) 
G[x(t),t] = — / — exp[ikx(T)} 2^ jt- 

(0.3) 

where, L is the chain length, a is the lattice constant, (3 
is the inverse temperature, v n are the fermionic Matsub- 
ara frequencies and = — Jcos(fc) is the electron dis- 
persion relation. The chemical potential has been pinned 
to the zero energy level. After summing over the frequen- 
cies in eq.(3) we observe that the average energy per site 
can be written as 

<H{ ^ ]> = V(x(r))+u(T)j(r) 

a r /a 

V(x(t)) = — J— / dkcos[kx(T)]cosh(ekT/fi)nF(ek) 
71 Jo 

a r /a 

j(r) = —a— / dkcos[kx(T)]cosh(ekT/h)riF(£k) 
* Jo 



(0.4) 

with N = L/a being the number of lattice sites and 
rip is the Fermi function. V{x{j)) is an effective term 
accounting for the r dependent electronic hopping while 
j(r) is the external source [20] current for the oscillator 
field u(t). Averaging the electrons over the ground state 
we neglect the fermion-fermion correlations [21] which 
lead to effective polaron-polaron interactions in non per- 
turbative analysis of the model [22]. This approximation 
however is not expected to affect substantially the fol- 
lowing calculations. 

Taking a large number of oscillators (ui(r),i = 1..N) 
as the bath for the quantum mechanical particle whose 
coordinate is x(t), one can write the general path integral 
at any temperature as: 



z(/j)ix(o) >= n/ d "'< t > / DxtT > ■ 

1 f hti (m 

(0.5) 

where, m is the electron mass, Mi is the mass of the 
i — th ionic group, u>i is the oscillator frequency, 7j is the 
coupling constant between external source term j(r) and 
i — th oscillator m. In our model, Mj = M and n = 1 
hence, integrating out the oscillators coordinates over the 
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paths Dui(r), imposing a closure condition on the paths 
and replacing t — > t/Ti, we get 



< 



x((3)\x(0) >= Y[Zi [ Dx{t)- 

i 



ft 2 A 1 



A ( X ( T V AM S 7^ sinh(?i^/3/2) 
jf drj(r) ^ rfr" cosh(^(|r - r"| - /3/2))j(r"), 

^ = 2sinh(^/3/2) ( °- 6) 

The quadratic source term A(x{t)) is the non local 
(in time) action for the particle having memory of the 
interactions with the oscillators field. Assuming periodic 
conditions x{t) = x(t + /3), the particle paths can be 
expanded in Fourier components 

oo 

x(t) = x Q + ^ 2^3?x„ cos(w„r) — 3x„ sin(u n T)^ 

n=l 

uj n = 2im/P (0.7) 
and, taking the following measure of integration 



u X (t)= r * x ° n 



< 2ttK /ttiKbT 



dSix n d^X n 

(0.8) 



a direct integration of eq. (6) can be carried out in order 
to derive the full partition function of the system versus 
temperature. 

Let's point out that, by mapping the electronic hop- 
ping motion onto the time scale, a continuum version of 
the interacting Hamiltonian (eq.(2)) has been de facto in- 
troduced. Unlike previous [23] approaches however, our 
path integral method is not constrained to the weak e-ph 
coupling regime and it can be applied to any range of 
physical parameters. 



II. COMPUTATIONAL METHOD AND RESULTS 



As a preliminar step we determine, for a given path and 
at a given temperature: i) the minimum number (Nk) 
of k— points in the Brillouin zone to accurately estimate 
the average interacting energy per lattice site and, ii) 
the minimum number (N T ) of points in the double time 
integration to get a numerically stable source action in 
eq.(6). We find: N k = 70 and N T = 300 at T = IK. 

Computation of eqs.(6)-(8) requires fixing two sets of 
input parameters. The first set contains the physical 
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quantities characterizing the system: the bare hopping 
integral J, the oscillator frequencies loi and the effective 
coupling x = a 2 h 2 /M (in units meV 3 ). The second set 
defines the paths for the particle motion which mainly 
contribute to the partition function through: the num- 
ber of pairs ($lx n , ^sx n ) in the Fourier expansion of eq.(7), 
the cutoff (A) on the integration range of the expansion 
coefficients in eq.(8) and the related number of points 
(N\) in the measure of integration which ensures numer- 
ical convergence. 

After introducing a dimensionless path 

x(t)/cl = x + ^2 (a n cos(w„t) + b n sin(w„r)^ , (0.9) 

71=1 

with: a n = 2$ix n /a and b n = — 23x„/a, the functional 
measure of eq.(8) can be rewritten for computational pur- 
poses as: 



r-A/a 



V2\2J ( 7 rh 2 /mK B T) N p+ 1 / 2 



2A/a 



da i 



2A/a 



2A/o 



dbx 



2 A/ a 



2A/o 



da Nv 



2A/a 



-A/a 
2A/a 



db 



2A/o 



(0.10) 



We take the lattice constant a = 1A. As a criterion 
to set the cutoff A on the integration range, we notice 
that the functional measure normalizes the kinetic term 
in eq.(6): 



<j> Dx{r)exp 



rn 
T 



Jo 



dTX 2 (r) 



1. 



(0.11) 



and this condition holds for any number of pairs N p 
truncating the Fourier expansion in eq.(9). Then, taking 
N p — 1 and using eq.(10), the left hand side transforms 
as: 



Dx(r)exp 



U 



A = 




mK B T 



drx 2 (r) 



u l2 
dyexp(-y 2 ) 



o 



(0.12) 



Using the series representation [24] 



/ dyexp(-y 2 ) = 



(— l) fe t/ 2fc+1 
A!(2fc + 1) ' 



(0.13) 



and taking [7=3 (the series converges with k max <~ 
100) eq.(13) yields 0.886207 which well approximates the 
Poisson integral value y/ir/2. 

Thus, the cutoff A can be expre ssed in terms of the 
thermal wavelength A as A ~ 3A/V27T 3 hence, it scales 
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versus temperature as A oc l/VT. This means physically 
that, at low temperatures, A is large since many paths are 
required to yield the correct normalization. For example, 
at T = IK, we get A <~ 284A. Numerical investigation of 
eq.(6) shows however that a much shorter cutoff suffices 
to guarantee convergence in the path integral, while the 
cutoff temperature dependence implied by eq.(12) holds 
also in the computation of the interacting partition func- 
tion. The thermodynamical results hereafter presented 
have been obtained by taking A <~ A/(10\/27r 3 ). Sum- 
ming over ~ 20/ VT points for each integration range 
and taking N p = 2, we are then evaluating the contribu- 
tion of (N\ + 1) 2N p +1 paths (the integer part of Na is ob- 
viously selected at any temperature). Thus, at T = IK 
we are considering <~ 4 • 10 6 different paths for the par- 
ticle motion while, at T = 100K the number of paths 
drops to 243. Low temperature calculations prove there- 
fore to be extremely time consuming. Note that larger 
N p in the path Fourier expansion would further increase 
the computing time without introducing any substantial 
improvement in the thermodynamical output of our cal- 
culation. 

Although the history of the SSH model is mainly re- 
lated to wide band polymers, we take here a narrow band 
system (J = lOOmeV) to be consistent with previous 
investigations [11] and with the caveat that electron- 
electron correlations may become relevant in narrow 
bands. Free energy and heat capacity have been first 
(Figures 1-3) computed up to room temperature assum- 
ing a bath of N = 10 low energy oscillators separated 
by 2meV: Jkvi = 2meV, huj w — 20meV. The low- 
est energy oscillator yields the largest contribution to the 
phonon partition function mainly in the low temperature 
regime while the wio oscillator essentially sets the phonon 
energy scale which determines the size of the e-ph cou- 
pling. A larger number N of oscillators in the aforegiven 
range would not significantly modify the calculation. 

In the discrete SSH model, the value a = 4a 2 / (ttkJ) ~ 
1, marks the crossover between weak and strong e-ph 
coupling, with n being the effective spring constant. In 
our continuum and semiclassical model the effective cou- 
pling is the above defined \. Although in principle, dis- 
crete and continuum models may feature non coincident 
crossover parameters, we assume that the relation be- 
tween a and J obtained by the discrete model crossover 
condition still holds in our model. Hence, at the crossover 
we get: %c ~ 7rJ7i 2 w 2 /64. This means that, in Figures 
1-3, the crossover value is set at Xc ~ 2000mey 3 . 

In Fig.l, the phonon free energy (F p h) is plotted ver- 
sus temperature together with the free energy due to the 
total action in eq.(6). F sou results from the competition 
between the free path action (kinetic term plus hopping 
potential in the exponential integrand) and the source 
action depending on the e-ph coupling. While the former 
enhances the free energy the latter becomes dominant at 
increasing temperatures thus reducing the total free en- 
ergy. As a main feature one notes that, by increasing 
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X, F sou gets a negative derivative at decreasing temper- 
atures. In the very weak coupling case (X = lOOmeV 3 ) 
F sou never intersects F p h whereas an intersection point 
shows up both at moderately weak (x = lOOOmeV 3 ) 
and at moderately strong (x = 3000meV 3 ) couplings. 
The intersection temperature decreases as expected by 
increasing the strength of \ but, at low T (-< 50K), the 
contribution of the free electron path action prevails: this 
feature is reflected on the low T behavior of the heat ca- 
pacity linear coefficient hereafter discussed. 

Fig. 2 shows the source term contributions to the heat 
capacity: the previously described summation over a 
large number of paths turns out to be essential to re- 
cover the correct thermodynamical behavior in the zero 
temperature limit. The small phonon contribution to the 
heat capacity is plotted in Fig. 3 to point out that the Du- 
long Petit value is achieved at T > 200K. Looking at the 
total heat capacity over temperature ratio, we find a pecu- 
liar low temperature upturn (also in the weak x regime) 
which can be mainly ascribed to the sizable effective hop- 
ping integral term V(x(t)). The e-ph coupling however 
determines the shape of the low T anomaly. 

The oscillators bath affects the electronic correlations 
both on the space and time scales. Analysis of the ther- 
mal correlation function shows infact that high energy 
phonons substantially reduce the electronic correlation 
length at low temperatures while, in presence of low en- 
ergy phonons, the electronic paths are correlated over the 
t scale. By increasing T, the role of the phonons on the 
correlation function is less pronounced. 

The effect of the oscillators bath on the thermodynam- 
ical properties is discussed in Figures 4-6 where the ten 
phonon energies are: huii = 22meV, ...,huJm = AOmeV. 
Accordingly the crossover is set at Xc ~ 8000mey 3 and 
three plots out of five lie in the strong e-ph coupling 
regime. As shown in Fig. 4, much larger x values (with re- 
spect to Fig.l) are required to get strongly decreasing free 
energies versus temperature while the x — 3000meV 3 
curve now hardly intersects the phonon free energy at 
room temperature. Fig. 5 shows the rapid growth of the 
source heat capacity versus temperature at strong cou- 
plings whereas the presence of the low T upturn in the 
total heat capacity over T ratio is confirmed in Fig. 6. 
Note that, due to the enhanced oscillators energies, the 
phonon heat capacity saturates at T <~ 400K . 

According to our integration method (eq.(10), at any 
temperature, a specific set of Fourier coefficients defines 
the ensamble of relevant particle paths over which the 
hopping potential V{x(t)) is evaluated. This ensamble 
is therefore T dependent. However, given a single set of 
path parameters one can monitor the V(x(t)) behavior 
versus T. It turns out that the hopping decreases (as ex- 
pected) by lowering T but its value remains appreciable 
also at low temperatures (■< 20K). Since the dr inte- 
gration range is larger at lower temperatures, the overall 
hopping potential contribution to the total action is rel- 
evant also at low T. It is precisely this property which is 
responsible for the anomalous upturn in the heat capac- 
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ity linear coefficient. Further investigation also reveals 
that the upturn persists both in the extremely narrow 
(J ~ lOmeV) and in the wide band (J ~ leV) regimes. 
Our computation method accounts for a variable range 
hopping on the r scale which corresponds physically to 
introduce some degree of disorder along the linear chain. 
This feature makes our model more general than the stan- 
dard SSH Hamiltonian (eq.(l)) with only real space near- 
est neighbor hops. While hopping type mechanisms have 
been suggested [25] to explain the striking conducting 
properties of doped polyacetylene at low temperatures 
we are not aware of any direct computation of specific 
heat in the SSH model. Since the latter quantity di- 
rectly probes the density of states and integrating over 
T the specific heat over T ratio one can have access to 
the experimental entropy, our method may provide a new 
approach to analyse the transition to a disordered state 
which indeed exists in polymers [26] . In this connection it 
is also worth noting that the low T upturn in the specific 
heat over T ratio is a peculiar property of glasses [27,28] 
in which tunneling states for atoms (or group of atoms) 
provide a non magnetic internal degree of freedom in the 
potential structure [29]. 



FIG. 1. Phonon and Source Term contributions to the 
free energy for three values of the effective coupling x (in 
units me V 3 ) and a narrow electron band. A bath of ten 
phonon oscillators is considered, the largest phonon energy 
being Tkdio = 20meV. 
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FIG. 2. Source Term contributions to the heat capacity for 
the same parameters as in Fig.l. 



FIG. 3. Total heat capacity over temperature for the same 
parameters as in Fig.l. The phonon heat capacity is also 
plotted. 
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FIG. 4. Phonon and Source Term contributions to the 
free energy for five values of the effective coupling \ (in 
units meV 3 ) and a narrow electron band. A bath of ten 
phonon oscillators has been taken, the largest phonon energy 
is hiOio = 40meV. 



FIG. 5. Source Term contributions to the heat capacity for 
the same parameters as in Fig. 4. 



10 



FIG. 6. Total heat capacity over temperature for the same 
parameters as in Fig. 4. The phonon heat capacity is also 
plotted. 

III. CONCLUSIONS 

Mapping the real space Su-Schrieffer-Heeger model 
onto the time space we have developed a semiclassical 
version of the interacting model Hamiltonian which can 
be attacked by path integrals methods. The acousti- 
cal phonons of the standard SSH model have been re- 
placed by a set of oscillators providing a bath for the 
electron interacting with the displacements field. Time 
retarded interactions are naturally introduced in the for- 
malism through the source action A{x{t)) which depends 
quadratically on the bare e-ph coupling strength a. Via 
calculation of the electronic motion path integral, the 
partition function can be derived in principle for any 
value of a thus avoiding those limitations on the e-ph cou- 
pling range which burden any perturbative method. Par- 
ticular attention has been paid to establish a reliable and 
general procedure which allows one to determine those in- 
put parameters intrinsic to the path integral formalism. 
It turns out that a large number of paths is required to 
carry out low temperature calculations which therefore 
become highly time consuming. The physical parameters 
have been specified to a narrow band system and the be- 
havior of some thermodynamical properties, free energy 
and heat capacity, has been analysed for some values of 
the effective coupling strength lying both in the weak and 
in the strong coupling regime. We find a peculiar upturn 
in the low temperature plots of the heat capacity over 
temperature ratio indicating that a glassy like behavior 
can arise in the linear chain as a consequence of a time 
dependent electronic hopping with variable range. 
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